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Abstract 

A complete discrete set of spherical single-particle wave functions for studies 
of weakly-bound many-body systems is proposed. The new basis is obtained 
by means of a local-scale point transformation of the spherical harmonic os- 
cillator wave functions. Unlike the harmonic oscillator states, the new wave 
functions decay exponentially at large distances. Using the new basis, charac- 
teristics of weakly-bound orbitals are analyzed and the ground state properties 
of some spherical doubly-magic nuclei are studied. The basis of the trans- 
formed harmonic oscillator is a significant improvement over the harmonic 
oscillator basis, especially in studies of exotic nuclei where the coupling to the 
particle continuum is important. 

PACS number(s): 21.10.Pc, 21.60.Cs, 21.60.Jz, 71.15.Mb 
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I. INTRODUCTION 



Understanding the nature of exotic nuclei with extreme isospin values is one of the 
most exciting challenges of current nuclear structure physics. Thanks to developments in 
radioactive ion beam instrumentation, we are in the process of exploring the very limits 
of nuclear binding, namely the regions of the periodic chart in the neighborhood of the 
particle drip lines. Several new structural features arise in the description of nuclei near the 
drip lines. Firstly, such nuclei exhibit very weak binding, which leads to extended spatial 
distributions. Secondly, the particle continuum plays a critical role in a description of the 
properties of these nuclei. Unlike more stable nuclei closer to the valley of beta stability, a 
proper theoretical description of weakly-bound systems (such as halo nuclei) requires a very 
careful treatment of the asymptotic part of the nucleonic density. 

The correct treatment of the particle continuum and near-threshold (e.g. halo) states 
poses a significant theoretical challenge. One possible way of tackling this problem is by 
means of the continuum shell model (CSM), in which the basis consists of both bound and 
unbound states, e.g., the eigenvectors of some finite shell-model potential An alter- 

native approach is to discretize the continuum by means of Sturmian function expansions 
or resonant state expansions. Sturmian functions form a discrete set of states which behave 
asymptotically as outgoing waves; they have been used as a basis in the solution of scat- 
tering equations, including various applications of the CSM @J|,§] . The Gamow (Berggren) 
states are the eigenstates of the time-independent Schrodinger equation with complex eigen- 



values 10-12]; they have been applied to many problems involving an unbound spectrum 



i~3Hl5|. Another possibility is to employ the canonical Hartree-Fock-Bogolyubov (HFB) 



basis of the independent-quasiparticle Hamiltonian ]IE|,r7||. The canonical states, i.e., the 
eigenstates of the one-body density matrix, form a complete localized basis with proper one- 
body asymptotics. Unfortunately, in order to find the canonical basis, the self-consistent 
one-body density matrix must be known, and this is not always possible. Finally, the par- 
ticle continuum can be discretized by placing the nucleus inside a very large box. Since 
the properties of the nucleus itself must not depend on the box size, one has to subtract 
the contribution from the free-gas states that are introduced ]T8|-pU[. The coordinate-space 
HFB [p!6| , |17| , pT|,p2| and relativistic Hartree-Bogolyubov (RHB) |23| methods are based on 
this principle. 

One of the most important tools of nuclear structure physics is the harmonic oscillator 
(HO) potential. The simple analytic structure of HO wave functions greatly simplifies shell- 
model studies f24|-p6||. Since the HO wave functions form a complete discrete set, they 
typically serve as the single-particle basis of choice in microscopic many-body calculations. 
Another useful feature of HO states is that they are all spatially localized, so that the 
resulting densities and currents are localized as well. The disadvantage of HO wave functions 
is that, because of their gaussian asymptotics, they cannot describe the falloff of nuclear 
density distributions near the nuclear surface and beyond. To get the falloff correct, it 
is essential to mix HO basis states. This is a serious problem when dealing with weakly- 
bound nuclei, where the density distributions fall off very slowly and the bases required are 
extremely large and thus very difficult to handle. Hence, the HO basis is not particularly 
useful when dealing with weakly-bound nuclei. 

It is desirable, therefore, to have an alternative to HO wave functions for use in nuclear 
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structure studies. Ideally, the new basis should preserve as many of the advantageous fea- 
tures of HO wave functions as possible, but, at the same time, it should have an improved 
asymptotic behavior. In this work, we consider a new spherical single-particle basis, ob- 
tained through the application of a general local-scaling point transformation to harmonic 
oscillator wave functions. This approach was originally developed by Petkov and Stoitsov, 
and is described in detail in Ref. [p7f| . In this study, we pay particular attention to the 
asymptotic form of the resulting transformed HO wave functions, developing the method 
so as to guarantee an appropriate exponential falloff. This makes the new basis especially 
useful in the description of weakly-bound nuclear systems. 

The paper is organized as follows. Section [TI| briefly reviews the method of Petkov and 
Stoitsov, and introduces the transformed HO basis. Stability tests of the new basis, when 
applied to weakly-bound orbitals, are contained in Sec. PTT| . In Sec. [TV], we use the new basis 
in variational calculations based on the density functional theory. Finally, Sec. [V| summarizes 
the principal conclusions of the work and spells out some issues for future consideration. 



II. TRANSFORMED OSCILLATOR BASIS 

A. Local-scale point transformations 

The key ingredient in the construction of our new basis is a coordinate transformation 
based on the local-scale transformation method [|27|. A local-scaling point transformation 



(LST) replaces the original coordinate r by a new coordinate r' = f(r) = r/(r). The new 
coordinate is in the same direction as r, but has a new magnitude r' = /(r), depending on 
a scalar function /(r) (called the LST function). It is assumed that /(r) is an increasing 
function of r and /(0)=0. The set of invertible transformations of this type forms a LST 
group. 

Given a model A-particle wave function ^(ri, r 2 , r^), the LST transforms it into a 
new wave function 



^/(ri,r 2 ,...,r A ) 



rf <9r j 



*(f(r 1 ),f(r a ),... l f(r x )). (1) 



Li=l 

Assuming that the model wave function is normalized to unity, 

<*!*) = 1, (2) 

the LST wave function ^/(ri, r 2 , r^) will also be normalized to unity, regardless of the 
choice of /(r). 

The local one-body density corresponding to an A-body wave function \1/ is 

p(r) = A J |^(r,r 2 ,...,r A )| 2 dr 2 ,...,dr4 . (3) 

It then follows from (TJ) that there exists a simple relation between the local density p/(r) 
associated with the LST function and the model local density p(r) corresponding to the 
model function ty: 
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(4) 



The relation (|J) is particularly useful when the density p/(r) is known (or at least approx- 
imately known). Given a model wave function, Eq. (|J) becomes a first-order nonlinear 
differential equation for the LST function /. For a spherically symmetric system, pf, p, and 
/ depend only on r=|r|, and Eq. (01) can be reduced to the nonlinear algebraic equation 



f(r) 



Pf(u)u du = / p(u)u du 



(5) 



which can be solved subject to the boundary condition /(0) = 0. Such an approach was 
pursued in a series of works [28-|30| based on the energy density functional method. A review 
of the density functional theory based on the LST can be found in the monograph ]31 . 



In the context of shell-model or mean-field applications, the particularly interesting case 
is when the model function is a Slater determinant, 



*(ri,r 2 , ...,r A ) 



det |^i(r 3 -)| 



3)\ ' 



(6) 



built from a complete set of model single-particle wave functions ipi(r). Due to the unitarity 
of transformation (|l|), the LST wave function retains the structure of a Slater determinant, 



r 2 , ...,r A ) 



but with new single-particle wave functions 



'A\ 



det |^i(rj)| 



(7) 



/ 2 (r) df(r) 

iy 2 (^f 



1/2 



UK*)) 



(8) 



These functions, in the following referred to as the LST basis, form a complete set of single- 
particle states. 



B. Transformed harmonic oscillator basis 



Having in mind the numerous advantages of the HO basis, we choose for the model 
single-particle wave functions the eigenfunctions of a HO potential. Since in this work we 
consider spherically symmetric systems only, the angular part of the single-particle wave 
function is not affected by the LST transformation. The radial HO wave functions R^f(r) 
are characterized by one external parameter, namely the oscillator length a osc = ^Ti/Mu. 
The LST basis associated with HO model wave functions will be referred to throughout the 
remainder of this paper as the transformed harmonic oscillator (THO) basis. The states in 
this basis are given by Eq. (H): 



^ H °(r) 



P{r) df(r) 



dr 



1/2 



*S?(/(r)). 



(9) 
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Up to this point, the LST function has not been defined. Ideally, we would like to 
parametrize f(r) in as simple manner as possible consistent with the requirement that it 
reproduces the generic features of the local density p/ at small and large distances r. Of 
course, the actual behavior of pf is affected by nuclear shell effects. Therefore, we concentrate 
on the average behavior of p/(r) only. 

In the nuclear interior, the average local density varies rather weakly with r. Due to 
the effects of Coulomb repulsion, the proton charge density is expected to exhibit a central 
depression p2|. Consequently, at low values of r, we assume the following ansatz: 



p(r) « p + cr . (10) 

Eq. (^) can then be cast into a fifth-order polynomial equation for the LST function f(r) at 
small r. To further simplify the problem, we assume that the model density p(r) is constant 
in the inner region. In this limit, the LST function assumes the simple form 

f(r) = r{a + br 2 ) 1/3 , (11) 

where a and b are parameters still to be determined. 

For large values of r, i.e., outside the nucleus, the average LST density should decay 
exponentially. Hence, 

p (r)^p(R)exp(^-) , (12) 



where R is a characteristic LST radius that is significantly greater than the nuclear radius 
and a is the LST diffuseness parameter. In contrast, the model density based on HO states 
exhibits a gaussian asymptotic behavior. Using Eq. (|^), we obtain an approximate expression 
for f{r) at large distances, 



/(r) = y ^ + ^ + d + d]nr + d x r , (13) 

where <i_2, cLi, do, d%, and d are parameters. Asymptotically, the linear term in the 
expansion ( |13"D takes over and f(r — > oo) ~ r 1 / 2 , as it should. 

The parameters {d-2, d-i, d\, do, d} can be determined by the requirement that f(r) and 
its first, second, third, and fourth derivatives are continuous at the point R. We are then 
left with only three independent parameters {a,b,R} in the LST function f(r), which can 
be chosen to optimize the THO basis for the physics problem of interest. In practical cal- 
culations, we assume that a, b > 0, thereby ensuring that f(r) is a monotonically increasing 
function of r. 



The ansatz (|ll]) -(|T3D for the LST function guarantees that all the THO states (§) are 
spatially localized and decay exponentially at large distances. The THO functions are con- 
tinuous up to their fourth derivatives. The HO length is absorbed into the coefficients a 
and b, and does not appear as an additional parameter. It should be stressed here that 
the parametrization of f(r) assumed in our work reflects our desire to keep it as simple 
and practical as possible. Although, as will be demonstrated in the following sections, this 
three-parameter form of the LST function performs very well in actual calculations, other 
choices are possible and could prove useful (or even essential) in other applications. 



5 



III. APPLICATION OF THE THO BASIS TO WEAKLY-BOUND STATES 



On the neutron-rich side of the valley of stability, there appear loosely-bound few-body 
systems called neutron halo nuclei (see Refs. [|33|-|3~6|1 for reviews). In these nuclei, weak neu- 
tron binding implies large spatial dimensions and the existence of the halo (i.e., a dramatic 
excess of neutrons at large distances). Theoretically, the weak binding and the correspond- 
ing proximity to the particle continuum, together with the need for explicit treatment of 
few-body dynamics, makes the subject of halos both extremely interesting and difficult. 

In this section, we apply the THO basis to weakly-bound single-particle states to assess 
its potential usefulness in the description of nuclei far from stability and, in particular, halo 
nuclei. We focus on spherical single-particle states and assume that they come from a finite 
square well (SQW) potential with radius Rq and depth — Vq. We concentrate on £=0 orbitals, 
since they are the best candidates for halos p7| , p8| . 

More specifically, we carry out a diagonalization of the SQW Hamiltonian within a 
truncated THO single-particle basis and compare the results that emerge with those of 
the exact SQW solutions. For the sake of comparison, we also carry out calculations in the 
analogous truncated HO basis. All calculations are carried out for a fixed radius i?o=7.11 fm 
of the SQW potential and under the assumption that the most weakly-bound £=0 state is 
the 3s. To simulate scenarios with varying degrees of binding of the 3s orbital, we vary 
the well depth. The parameters defining the basis, i.e., the oscillator length in the HO 
variant and the {a, b, R} parameters of the THO basis, have been chosen so as to minimize 
the single-particle energy of the 3s halo state. For a given calculation in either the HO or 
THO basis, truncation is defined to include all single-particle states belonging to iV<iV max 
oscillator shells. 

The ability of the HO and THO basis expansions to reproduce single-particle energies 
is illustrated in Fig. [I], which shows the deviation between approximate and exact energies 
of the Is, 2s, and 3s states as a function of iV max . The energy of the 3s state was assumed 
in these calculations to be very low: —200 keV (Fig. |I]a) and — 40keV (Fig. |l]b). As can be 
seen from the figure, the THO basis offers a systematic improvement over the traditional 
HO expansion. As an example, iV max =20 THO shells are sufficient to reproduce the energy 
of the 3s halo state with an accuracy of 50-60 keV, whereas iV max ^30 is required with the 
HO basis. For the well-bound Is and 2s states, there is also an improvement when using 
the THO basis. (The low-iV max fluctuations seen in the Is and 2s curves reflect the fact 
that the basis was optimized to the energy of the 3s halo state, so that the resulting radial 
asymptotics is not appropriate for the more deeply-bound orbitals.) 

Figures [§|4] compare the exact and approximate wave functions. It is gratifying to see 
that even with a relatively low number of THO shells the structure of the 3s orbital (Figs. 
H and |3|) is well reproduced out to 15-18 fm, and the agreement with the exact eigenvector 
becomes excellent for iV max =30. Again, as for the single-particle energies, one needs at least 
30 HO shells to obtain results of comparable quality out to 15-18 fm. For well-bound states, 
such as the Is state at -22.4 MeV shown in Fig. |], both expansions work equally well: 20 
HO or THO shells are sufficient to reproduce the exact result in the physically interesting 
region. 

It is instructive to discuss why so many HO shells are required to reproduce a wave 
function out to a very large distance. In the HO approximation, the classical radius of an 
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orbit with principal quantum number N and orbital quantum number £=Q is given by 



r d « a osc V2N . (14) 

From this we see that the radial information contained in HO wave functions varies quite 
slowly with the principal quantum number. To build up the large-r dependence in a radial 
wave function, we must include HO states of very high principal quantum numbers. In 
contrast, the THO basis has no such restriction on the radial content of its wave functions 
and convergence can be achieved much more rapidly. 

Another rather extreme example is shown in Fig. |5|. Here, we consider the single-particle 
wave function of a Is state at e=— 14keV. The SQW wave function has a very large spatial 
extension. In this case, instead of performing the full basis expansion, the optimized single- 
particle Is HO and THO wave functions were obtained by maximizing their respective 
overlaps with the exact solution. (In this way one is testing the ability of unperturbed basis 
states to reproduce the exact wave functions of weakly-bound states.) The corresponding 
wave functions are plotted in Fig. [5[ The THO wave function has a squared overlap with 
the exact SQW solution of | (SQW|THO) | 2 =0.994, whereas the HO wave function has a 
much lower squared overlap of | (SQW|HO) | 2 =0.877. Clearly, the THO basis with the LST 
functions (p]) -([13|) is much better able to reproduce the tail of a halo wave function than 
the HO basis. 



IV. VARIATIONAL CALCULATIONS WITH SKYRME FORCES 

In this section, we test whether Slater determinants built up in terms of the single- 
particle HO and THO wave functions are able to reproduce the results of full self-consistent 
HF calculations. We follow the energy density functional approach |28|j3T|j40[| , whereby the 



total HF energy, taken as an expectation value of the nuclear Hamiltonian over a trial Slater 
determinant, involves a sum of the Skyrme and Coulomb energies. As usual, the Skyrme 
energy density is expressed in terms of local nucleon densities, kinetic energy densities, 
and spin-orbit densities, all defined in terms of the variational single-particle states. The 
Coulomb energy density, which depends on the local proton density, contains both direct 
and exchange terms, the latter taken in the Slater approximation. The effective interaction 



used in these calculations was the Skyrme force SkP fl6 |. 

In the HO variational analysis, the energy functional of the SkP Hamiltonian was mini- 
mized with respect to two parameters, the proton and neutron harmonic-oscillator lengths. 
In the THO analysis, the energy minimization involved the six parameters that define the 
LST functions for neutrons {a n ,b n , R n } and protons {a p ,b p , R p }. The calculations were 
carried out for three spherical doubly-magic nuclei, 16 0, 40 Ca, and 208 Pb. 

Table | contains the results for ground-state binding energies and proton and neutron 
rms radii. As expected, the binding energies based on THO Slater determinants are lower 
than those based on the HO. The HF binding energies are of course lower than those based 
on either the HO or THO Slater determinants because of the self-consistent nature of HF 
calculations (the HF state is the optimal Slater determinant for a given Hamiltonian). In all 
cases considered, the THO binding energies are within 1.1% of the HF results, while the HO 
results deviate up to 5% for heavy nuclei. For neutron and proton rms radii, the THO results 
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agree very well with the self-consistent HF values, whereas the HO results systematically 
overestimate them. 

Table [TI] compares the neutron single-particle energies that resulted from the restricted 
HO and THO variational calculations with those from the self-consistent HF calculations. 
The restricted calculations give a rather good approximation to the self-consistent single- 
particle energies. In general, the best agreement is obtained for the high-£ states where the 
tail of the wave function plays little role. The largest deviations are seen for s and p states. 

Finally, the neutron and proton density distributions resulting from the various cal- 
culations are compared in Fig. |5|. There is excellent agreement between the THO and 
self-consistent HF densities in the surface region. In contrast, the incorrect asymptotic be- 
havior of HO wave functions leads to significant deviations |41 |. The THO densities differ 
most significantly from the HF results at small values of r. However, since the densities are 
weighted by r 2 in calculations of expectation values, the main contribution to global nuclear 
characteristics such as energies and rms radii comes from the surface region. 



V. SUMMARY AND CONCLUSIONS 

In this paper, we have explored a new class of single-particle basis states obtained by 
a local-scale point transformation (LST) of harmonic oscillator states. We focussed special 
attention on the asymptotic properties of these states (called THO states), to see whether 
they might be useful in the description of weakly-bound nuclear systems, including those 
with a halo structure. 

Following a comprehensive summary of the LST formalism and its use in building the 
THO basis, we discussed two applications of this new basis. The first concerned the de- 
scription of sub-threshold (halo) states. We showed that the THO basis is greatly superior 
to the usual HO basis in reproducing the properties of such weakly-bound states. We then 
discussed the use of this basis in restricted HF calculations. Once again, the clear superiority 
of the THO basis to the ordinary HO basis was demonstrated, even for normal well-bound 
nuclei. Most importantly, the optimal THO basis that emerged from the restricted HF 
calculations provided an excellent reproduction of nuclear surface properties. 

The analysis presented in this paper should be viewed as a starting point for future 
investigations. One potentially interesting application of the THO basis is in the context of 
nuclear shell-model studies. The fact that the THO basis expansion technique and the THO 
variational procedure accurately reproduce nuclear properties in the surface region suggests 
that the new basis can be very useful when studying those nuclear properties that depend 
on the asymptotic behavior of the radial form-factor, and also for microscopic calculations of 



the effective interactions for weakly bound systems [Q. Some of the simplicity inherent in 
the use of HO wave functions in such studies will be lost, however. Calculations of G- matrix 
elements, for example, benefit greatly from the ease of transforming harmonic oscillator 
product wave functions into relative and center-of-mass coordinates, and this would be lost 
in the THO basis. On the other hand, it should be still much easier to carry out calculations 
in this basis than in a basis generated by numerical solution of the single-particle Schrodinger 
equation. 

Another interesting avenue for future exploration concerns the use of the THO basis in 
HFB or RHB calculations of weakly-bound nuclei where other traditional methods (e.g., 
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HF+BCS) cannot be used [JT6l,17 . The fact that with a rather modest number of basis 
THO states one can reproduce the properties of spatially extended states suggests that 
the traditional method of solving the HFB equations, based on basis expansion, can be 
revitalized by using the THO basis. This way of solving the HFB problem (both in the 
spherical and deformed cases) can be an interesting alternative to algorithms based on 
coordinate-space methods |2l|-p3|. 

Finally, the THO states that arise variationally can be used as approximations to the 
self-consistent canonical states of HFB. Recently, the variational method based on the energy 
density formalism was generalized to the HFB case [1"3|. In that work, the authors employed 
the basis proposed by Ginocchio in the description of semi-magic nuclei. Despite the 
successes achieved in those calculations and in earlier restricted HF calculations for doubly- 
magic nuclei [45|| , there are several features of the Ginocchio potential that make its more 
detailed use problematical. First, since the centrifugal term in the Ginocchio potential is 
not treated properly, the resulting wave functions with £>0 do not have proper asymptotic 
behavior. Second, since the Ginocchio potential is finite, its bound eigenvectors do not form 
a closed set, and the use of the continuum wave functions is necessary (see recent Ref. [p5fl ). 
In this context, the discrete THO basis would appear to be more useful. 
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TABLES 



TABLE I. Total binding energies and neutron and proton rms radii obtained using HO and 
THO Slater determinants and from self-consistent HF calculations for the spherical doubly-magic 
nuclei 16 0, 40 Ca, and 208 Pb. 







16q 






40 Ca 






208p b 






HO 


THO 


HF 


HO 


THO 


HF 


HO 


THO 


HF 


E (MeV) 


-125.9 


-126.2 


-127.6 


-339.1 


-340.9 


-343.3 


-1565.5 


-1621.8 


-1636.5 


(rp)rms (fm) 


2.73 


2.73 


2.73 


3.46 


3.45 


3.45 


5.49 


5.47 


5.47 


(rn)rms (fm) 


2.71 


2.71 


2.70 


3.42 


3.40 


3.40 


5.64 


5.61 


5.61 



TABLE II. Spherical single-neutron energies (in MeV) obtained using HO and THO Slater 
determinants and from self-consistent HF calculations for the spherical doubly-magic nuclei 16 O, 
40 Ca, and 208 Pb. 



nlj 


HO 


16Q 

THO 


HF 


HO 


40 Ca 
THO 


HF 


HO 


208p b 

THO 


HF 


lSl/2 


-29.0 


-29.2 


-29.0 


-36.4 


-37.1 


-37.4 


-37.1 


-39.0 


-40.0 


1^3/2 


-17.9 


-17.9 


-18.3 


-29.0 


-29.1 


-28.9 


-34.8 


-35.6 


-36.4 


lfl/2 


-13.5 


-13.6 


-13.6 


-25.5 


-25.5 


-25.8 


-33.3 


-35.1 


-35.8 


14/2 








-19.6 


-19.6 


-19.7 


-31.3 


-31.6 


-31.9 


2si/ 2 








-14.9 


-15.5 


-16.4 


-28.7 


-29.8 


-29.2 


14/2 








-14.5 


-14.5 


-14.4 


-29.0 


-30.4 


-30.6 


L/7/2 














-26.7 


-26.7 


-26.9 


2^3/2 














-23.2 


-23.8 


-23.3 


1/5/2 














-23.5 


-24.5 


-24.5 


2pi/ 2 














-21.7 


-22.7 


-22.4 


159/2 














-21.1 


-21.1 


-21.3 


24/2 














-17.1 


-17.5 


-17.2 


lfiV/2 














-17.0 


-17.6 


-17.6 


3«l/2 














-15.0 


-15.8 


-15.5 


2 4/2 














-14.6 


-15.5 


-15.6 


1^11/2 














-14.8 


-15.0 


-15.3 


2/7/2 














-11.0 


-11.3 


-10.9 


lhg/2 














-9.7 


-10.0 


-10.1 


1*13/2 














-8.0 


-8.5 


-8.8 


3^3/2 














-9.0 


-9.5 


-8.9 


2/5/2 














-7.8 


-8.5 


-8.6 


3pi/ 2 














-7.7 


-8.4 


-8.1 
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FIGURES 



FIG. 1. Deviation between approximate energies based on the basis expansion method and 
exact energies of the Is, 2s, and 3s states of the SQW Hamiltonian as a function of the number of 
HO quanta included in the basis, N max . The solid lines are the results obtained when expanding 
in a HO basis; the dotted lines are the results obtained using the THO basis. The radius of the 
well was taken to be Rq=7.11 fm. The different binding energies of the 3s halo orbital, —200 keV 
(a) and — 40keV (b), were achieved by changing the well depth. 

FIG. 2. Exact (dots) and approximate (lines) wave functions of the 3s halo state at e=— 200 keV 
as a function of r. The approximate wave functions were obtained by means of the basis expansion 
method (HO, left panel; THO, right panel) in iV max =12, 20, 30, and 40 HO (THO) shells. 

FIG. 3. Same as in Fig. [2] but for the wave functions of the 3s halo state at e=— 40 keV and for 
iV max =12, 20, 30, and 60 HO (THO) shells. 

FIG. 4. Same as in Fig. | but for the wave functions of the well-bound Is state at e=— 22.4MeV. 

FIG. 5. The wave function of a Is SQW eigenstate at e=-14keV (SQW, solid line). The 
HO (dotted line) and THO (dashed line) Is wave functions were determined by maximizing their 
overlap with the SQW state. The asymptotic behavior (in logarithmic scale) is shown in the insert. 

FIG. 6. Local one body densities for 16 0, 40 Ca, and 208 Pb calculated in the self-consistent 
HF approach (solid line) compared with those obtained in restricted variational calculations using 
as trial wave functions Slater determinants built from HO (dotted line) and THO (dashed line) 
single-particle orbitals. All calculations utilized the same SkP effective interation. 
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